Paths of fluctuation induced switching 
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We demonstrate that the paths followed by a system in fluctuation-activated switching form a 
narrow tube in phase space. A theory of the path distribution is developed and its direct measure- 
ment is performed in a micromechanical oscillator. The experimental and theoretical results are in 
excellent agreement, with no adjustable parameters. We also demonstrate the lack of time-reversal 
symmetry in switching of systems far from thermal equilibrium. 
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Fluctuation-activated switching between coexisting 
stable states is at the root of diverse phenomena, from 
switching in nanomagnets [l| and Josephson junctions 
Q to chemical reactions and to protein folding jj] • A de- 
tailed theory of switching rates was first developed by 
Kramers [H for systems close to thermal equilibrium. 
Here, the switching rate is determined by the free en- 
ergy barrier between the states. In recent years much 
attention has been given to switching in systems far from 
thermal equilibrium, like electrons Q and atoms 0, 
modulated traps and rf-driven Josephson junctions 
and nano- and micromechanical resonators (lol . [ill . IK 
Nonequilibrium systems generally lack detailed balance, 
and the switching rates may not be found by a simple 
extension of the Kramers approach. 

The analysis of switching in both nonequilibrium 
systems and complex equilibrium systems, including 
biomoleculcs, relies on the idea that, even though the 
system motion is random, the switching paths form nar- 
row tubes in phase space. The tube is centered at the 
most probable switching path (MPSP). For low fluctua- 
tion intensity, the MPSP is obtained from a variational 
problem, which also determines the switching activation 
barrier 0, [lj, H [lj, [U H [H H3, Ell, SI ■ Despite 



its fundamental role, the concept of the narrow tube of 
switching paths has not been tested experimentally nor 
has this tube been characterized quantitatively. Prior ex- 
perimental studies [23|] and simulations 24 , 25 1 focused 
on the distribution in space and time of fluctuational 
paths to a certain space-time point [2tSj . The methods 



231 . 24 , 25 , |26( do not apply to the switching paths dis- 
tribution for systems with more than one dynamical vari- 
able because of the motion slowing down near the stable 
states and the related loss of path synchronization. 

In this paper, we present a theory of the switching 
paths distribution and report its direct observation in a 
well-characterized system, a micromechanical oscillator 
driven into parametric resonance. The observed distribu- 
tion is shown in Fig. [IJl. Our experimental and theoreti- 
cal results are in excellent agreement, with no adjustable 
parameters. In addition to proving the concept of the 
narrow tube of switching paths and putting it on a quan- 
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FIG. 1: (a) Switching probability distribution in a paramet- 
rically driven microelectromechanical oscillator. The distri- 
bution pi2(X, Y) is measured for switching out of state A\ 
into state Ai. (b) The peak locations of the distribution 
are plotted as black circles and the theoretical most probable 
switching path is indicated by the red line. All trajectories 
originate from within the green circle in the vicinity of A\ and 
later arrive at the green circle around A2- The radii of the 
green circles give the typical fluctuation amplitude Id- The 
portion of the distribution outside the blue lines is omitted, 
(c) Scanning electron micrograph of the micromechanical tor- 
sional oscillator, (d) Cross-sectional schematic. 



titative basis, the results provide the first demonstration 
of the lack of time-reversal symmetry in switching of sys- 
tems far from thermal equilibrium. The results also open 
the possibility of efficient control of the switching proba- 
bility based on the measured narrow path distribution. 

We consider a bistable system with several dynami- 
cal variables q = (qi, ...,qjy). The stable states A\ and 
A 2 are located at and (\a 2 respectively. For low 
fluctuation intensity, the physical picture of switching is 
as follows. The system prepared initially near state A±, 
for example, would fluctuate about it with small typi- 
cal amplitude Id. Ultimately, over time of order of the 
reciprocal switching rate WjT, 1 that largely exceeds the 
relaxation time t r , the system will make a transition to 
state A2 and will then fluctuate about it. In the transi- 
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tion the system first moves from the vicinity of qA 1 to 
the vicinity of the saddle point qs as a result of a large 
and rare fluctuation. Its trajectory is expected to be 
close to the one for which the probability of the appro- 
priate fluctuation is maximal. From the vicinity of qs 
the system moves to state A2 close to the deterministic 
fluctuation-free "downhill" trajectory. These two trajec- 
tories comprise the MPSP. 

We characterize the switching paths distribution by 
the conditional probability density for the system to pass 
through a point q on its way from Ai to A2, 

Pi2(q,£)=/ dq f p(qf,t f ;q,t\q ,to). (1) 

Here, the integrand is the probability density for the sys- 
tem to be at points q/ and q at times t/ and t, respec- 
tively, given that it was at qo at time to- The point qo 
lies within distance ~ Id of q^ and is otherwise arbi- 
trary. Integration with respect to q/ goes over the range 
f^2 of small fluctuations about q^ 2 ■ 

For Ww 1 , W21 ~>tf—t, t-t > t r , far from the stable 
states the distribution |T]) is independent of qo, to, tf, 
and t. It peaks on the MPSP. The peak is Gaussian 
transverse to the MPSP. For |q — q^^l, |q — qs| 3> Id 

Pi2(q,t)-^ 12 ^ 1 (e||)Z- 1 exp^-^ ± Q^j , (2) 

where £|| and £j_ are coordinates along and transverse 
to the MPSP, and v(^) is the velocity along the MPSP. 

The matrix elements of matrix Q ~ are oc l~^ 2 , and 

Z= [(2^) Ar - 1 /detQ] 1 /2. 

We derive Eq. @ for a system described by the 
Langevin equation of motion 

q = K(q) + f(t), (f n (t)f m (t')) = 2DS nm S(t - t'). (3) 

Here, f (t) is white Gaussian noise. The noise intensity D 
is assumed small. The dependence of the switching rates 
Wjj on D is given by the activation law, logWjj oc D^ 1 
0. The length l D = (DU) 1 / 2 where the relaxation 
time t r is given by the reciprocal minimal (in absolute 
value) real part of the eigenvalues of the matrix dK m /dq n 
calculated at q^i^: where K(q) = 0. 

The integrand in Eq. (1) can be written as a prod- 
uct of transition probabilities, p(q/, i/; q, t|qo, to) — 
p(q/, t/|q, t)p(q, t|qo, to)- In the time range of interest 
they satisfy the stationary backward and forward Fokker- 
Planck equations, respectively, 

[Kd^ + Ddl] p(q f ,t f \q,t) = 0, 

[-0 q K + £>^]Kq,t|qo,to) = O. (4) 

We start with the case where the observation point 
q lies in the region of attraction of the initially empty 
state A 2 . Here, p(q/,t/|q, t) sa p 2 (q.f), where p 2 (q/) is 



the normalized quasistationary probability distribution 
in the state A 2 in the neglect of interstate switching. 
Then, from Eq. (1) Pi 2 (q,t) = p(q, t|q , t ). 

The function p(q, t|qo, to) describes a stationary prob- 
ability current from the vicinity of qs to the attractor 
q^4 2 . It is concentrated on the noise-free "downhill" tra- 
jectory q = K coming out of the saddle, which thus gives 
the corresponding section of the MPSP. Diffusion leads 
to broadening of the distribution p(q, t|q ,tp) [27j | . The 
total current gives the switching rate Wi 2 [4(. 

On the downhill trajectory, ^ points along the vector 
K and the velocity is u(£||) = K(£n,£j_ = 0). For small 
£j_ < Id the solution of Eq. (j4|) has the form (2), with 
matrix Q given by equation v Q + k' i Q + Qk + 2Q 2 D = 

0, where k^ v = dK^/dt;± v with the derivatives evaluated 
for £j_ = 0; the subscripts \i, v = 1, . . . , N — 1 enumerate 
the components of £j_ and the transverse components of 
K in the co-moving frame. 

The case where the observation point q lies in the basin 
of attraction of the initially occupied state A\ is more 
complicated. We seek the solution of Eqs. for small 
D in the eikonal form, p(q, t|q , to) = exp[— Sp(q)/D], 
P(q/!*/|q>*) = exp[5 , B (q)/ J D]/o 2 (q/), with 5 f ,b = 
sf B + DS^ B + .... To the lowest order in D we have 

H (q, d q 4 0) B ) = 0, iJ(q, p) = p 2 + P K(q). (5) 

Equation (|5|) has a form of a Hamilton-Jacobi equa- 
tion. The function ^^(q) is the least action for reach- 
ing q starting from q^ []]|; S^^q) is the least ac- 
tion for reaching qs from q. The MPSP is given by 
the heteroclinic Hamiltonian trajectory from q^ to qs. 
Clearly, on the MPSP S B 0) (q) - sf\q) is maximal, 

1. e., <9 q 4 0) = dqSg\ From Eq. ©, the MPSP di- 
rection £|| and the velocity are given by the expression 

29 q S F 0) B(q)+K(q)= v(^)i {] . 

Matrix Q in Eq. ^ can be found from Eq. ©, Q = 

Qf - Qb, (Qf,bW = D- 1 d 2 S^JdU^ ±v with the 
derivatives calculated on the MPSP, i.e., for = 0. The 
prefactor in p\2 (q, t) is determined by — S B . Central 
for obtaining it is the relation vc^ Q + 2D ^Qf — Q^j + 

k^Q + Qk = that follows from Eq. (HJ). 

An important property of the distribution is prob- 
ability current conservation, J d^±v(^)pi2(q,t) = W\ 2 
independent of the position on the path. This result 
could be expected from the physical picture of switching. 
Because the velocity v(£||) decreases near the stationary 
states qAi, 2 , qs, the distribution increases there. Close 
to the saddle point Eq. @ should be modified. For sys- 
tems with detailed balance one should replace v by Qd /tr 
with £ ~ 1. 

The switching probability distribution ((T|) is measured 
for a high-Q micro-electro-mechanical torsional oscillator 
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(Q = 9966) driven into parametric resonance [28|. As 
shown in Fig. [TJ:, the device consists of a 3.5-/xm thick, 
200-//TO 2 heavily doped polysilicon plate suspended by 
two thin torsional rods on opposite edges. In addition 
to the restoring torque of the torsional springs, the top 
plate is also subjected to an electrostatic torque when a 
voltage is applied to the left electrode in Fig. [TJi. The 
voltage is modulated at a frequency lu/2tt (44346.800 Hz) 
close to twice the natural oscillation frequency (22178.578 
Hz) of the plate. The orientation angle of the plate 9(t) 
is detected capacitively through the right electrode in 
Fig. [TH. 

When the modulation is sufficiently strong, the plate 
oscillates at half the modulation frequency as a result of 
parametric resonance. Since the modulation is invariant 
upon a shift in time by its period, there exist two stable 
oscillation states that have the same amplitude but differ 
in phase by n. The dynamics is well described by the 
rotating wave approximation [28| . It is characterized by 
two dynamical variables X and Y (the quadratures), with 
9(t) = X(t) cos(u>t/2)+Y(t) sin(arf/2). A lockin amplifier 
is used to record X and Y every 2 ms. In Fig. [TJo, we 
show the two stable oscillation states, Ai and A-i- They 
are located symmetrically about the origin in the (X,Y)- 
plane. The measurement uncertainty is ~ 80 /irad, about 
0.6% of the full scale in Fig. [TJs. All measurements were 
performed at 77 K and < 10~ 6 torr. 

When white noise is added to the excitation voltage, 
the system can occasionally overcome the activation bar- 
rier and switch from one stable state to the other. The 
noise intensity is chosen to ensure that the mean resi- 
dence time in each state (~ 10 s) is much larger than 
the relaxation time (~ 1 s) of the system. The os- 
cillator dynamics in the rotating frame is described by 
Eq. © with <7i = CX, q 2 = CY and K = -C -1 q + eVg, 
where g = g 4 /4 - (1 - m)?i/ 2 + (1 + m)*?!/ 2 and ^ is tlie 
permutation tensor; from the harmonic and parametric 
resonances of the oscillator we found C = 176.349, £ = 
4.968, fi= 0.9367. 

Transitions are identified when the oscillator begins 
in the vicinity of the stable state A\ (within the left 
green circle in Fig. [T}j) and subsequently arrives at state 
A2 (within the right green circle). Figure [TJi shows 
the switching probability distribution derived from more 
than 6500 transitions. While in each transition the sys- 
tem follows a different trajectory, the trajectories clearly 
lie within a narrow tube. The maximum of the distri- 
bution gives the MPSP. In Fig. [TJo, the location of this 
maximum is plotted on top of the MPSP obtained from 
theory. The MPSP emerges clockwise from A\ and spi- 
rals toward the saddle point at the origin. Upon exit 
from the saddle point, it makes an angle and continues to 
spiral clockwise toward A2 . There is excellent agreement 
between the measured peak and the MPSP obtained from 
theory, with no adjustable parameters. 

Close to the stable states the peaks of the distribu- 



tion at successive turns of the MPSP overlap, preventing 
the accurate determination of the MPSP. The plot in 
Fig. []}L,b has excluded the portions of trajectories prior 
to escaping from the initial state A-y and upon arriving 
at the final state A%, which are bound by the two blue 
lines. Such cutoff also eliminates the large peaks of the 
distribution centered at A\ and A2, which arise because 
the oscillator spends most of its time fluctuating about 
Ai and A2. These peaks are not relevant to switching 
dynamics. 

The switching probability distribution in our multi- 
variable system displays important generic features. Fig- 
ure [5Ji shows the distribution cross-section along the pur- 
ple line transverse to the MPSP in Fig. [TJj. It is well- 
fitted by a Gaussian. Gaussian distributions with dif- 
ferent height and area are observed also in other cross- 
sections except close to the saddle point. Figure [2}j plots 
the area under the Gaussian distribution versus the re- 
ciprocal measured velocity on the MPSP, for different 
cross-sections. The linear dependence agrees with Eq. ([2]) 
and indicates that the probability current from the ini- 
tially occupied attractor to the empty one is constant. 
Near the saddle point the velocity is significantly dimin- 
ished and the motion becomes diffusive, leading to strong 
broadening and increase in height of the distribution. 
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FIG. 2: (a) The cross-section of the switching probability 
distribution along the purple line in Fig. [Tb transverse to the 
MPSP. The solid line is a Gaussian fit. (b) Velocity on the 
MPSP vs. inverse area under cross-sections of the switching 
probability distribution. The solid line is a linear fit forced 
through the origin. 

We find that the probability current concentrates 
within a narrow tube deep into the basins of attraction of 
A\ and A2 . In the basin of attraction to A2 but not too 
close to A2, much of the probability distribution carries 
the switching current. However, the probability distribu- 
tion deep inside the basin of attraction of A\ is largely 
associated with fluctuations about A\ that do not lead 
to switching. The part of the distribution responsible for 
the switching current is an exponentially small fraction 
of the total distribution. Nevertheless our formulation 
makes it possible to single out and directly observe this 
fraction. 

Another generic feature of the observed distribution is 
characteristic of systems far from thermal equilibrium. 
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For equilibrium systems, the most probable fluctuational 
path "uphill" from an attractor to the saddle point is the 
time reversal of the fluctuation-free "downhill" path from 
the saddle point back to the attractor. More precisely, it 
corresponds to the change of the sign of dissipation term 
in the equation of motion, cf. Ref. [2^. Our paramet- 
ric oscillator is far from thermal equilibrium and lacks 
detailed balance. Upon sign reversal of dissipation, even 
the attractors are shifted away from the original locations 
(from A\ to A\ in Fig. [3|). Our data show that the uphill 
section of the MPSP is distinct from both the dissipation- 
reversed path and the downhill noise-free path from the 
saddle to the stable state. 



"D 

> 

















A S ^\ 

\ /\ ' 




/ H - a / 





" J 1 ' 







-12 



-8 -4 
X (mrad) 



FIG. 3: (Color online) Comparison of the MPSP, the time- 
reversed path and the deterministic downhill path. The sec- 
tion of the most probable switching path from A\ to S is 
shown as a thick solid line. Upon changing the sign of dis- 
sipation, the attractor is shifted to a new location A[ and 
becomes an unstable state. The fluctuation-free path with 
reversed dissipation from A[ to S is shown as the thin solid 
line. The dashed line represents the deterministic downhill 
path from S to Ai. 

The observation of the MPSP reported here provides, 
in some respects, an experimental basis for the broadly 
used concept of a reaction coordinate, which can be as- 
sociated with the coordinate along the MPSP. An ad- 
vantageous feature of our approach is that it does not 
rely on a specific model of the fluctuating system. The 
only characteristics used are the positions of the stable 
states in phase space, which are usually accessible to di- 
rect measurement. The approach applies to systems far 
from thermal equilibrium as well as to equilibrium sys- 
tems. Measuring the switching trajectories can help to 
determine the model globally, far from the stable states. 
It can also provide an efficient way of controlling the 
switching rates by affecting the system locally on the 
most probable switching path. 

In conclusion, we propose a theory and report the first 
observation of the distribution of paths followed by a sys- 
tem in switching between stable states. We show that, 
in phase space, switching paths form a narrow tube cen- 



tered at the MPSP. The observed features of the paths 
distribution and the position of the MPSP for a well- 
characterized system are in full agreement with the the- 
ory. 

This research was supported in part by NSF DMR- 
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[1] 
[3 



[4 
[5 

[6] 
[7] 
[« 
[9 
[10 

[11 

[12 
[13 

[14 

[15 
[16 

[17 
[18 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[28 

[29] 



W. Wernsdorfer et ai, Phys. Rev. B 55, 11552 (1997). 
T. Fulton and L. Dunkelberger, Phys. Rev. B 9, 4760 
(1974). 

D. Wales, Energy Landscapes: Applications to Clusters, 
Biomolecules and Glasses (Cambridge University Press, 
2003). 

H. Kramers, Physica (Utrecht) 7, 284 (1940). 

L. J. Lapidus, D. Enzer, and G. Gabrielse, Phys. Rev. 
Lett. 83, 899 (1999). 

R. Gommers et al, Phys. Rev. Lett. 94, 143001 (2005). 
K. Kim et al, Phys. Rev. Lett. 96, 150601 (2006). 

I. Siddiqi et ai, Phys. Rev. Lett. 93, 207002 (2004). 
A. Lupascu et al, Nature Physics 3, 119 (2007). 

J. S. Aldridge and A. N. Cleland, Phys. Rev. Lett. 94, 
156403 (2005). 

C. Stambaugh and H. B. Chan, Phys. Rev. B 73, 172302 
(2006). 

R. Almog et ai, Appl. Phys. Lett. 90, 013508 (2007). 
M. I. Freidlin and A. D. Wentzell, Random Perturba- 
tions of Dynamical Systems (Springer- Verlag, New York, 
1998), 2nd ed. 

M. I. Dykman and M. A. Krivoglaz, Zh. Eksp. Teor. Fiz. 
77, 60 (1979). 

R. Graham and T. Tel, Phys. Rev. Lett. 52, 9 (1984). 
A. J. Bray and A. J. McKane, Phys. Rev. Lett. 62, 493 
(1989). 

R. S. Maier and D. L. Stein, Phys. Rev. Lett. 71, 1783 
(1993). 

S. Kraut and C. Grebogi, Phys. Rev. Lett. 93, 250603 
(2004). 

M. I. Dykman, B. Golding, and D. Ryvkine, Phys. Rev. 

Lett. 92, 080602 (2004). 

M. Bier, Phys. Rev. E 71, 011108 (2005). 

O. A. Tretiakov and K. A. Matveev, Phys. Rev. B 71, 

165326 (2005). 

V. Elgart and A. Kamenev, Phys. Rev. E 74, 041101 
(2006). 

J. Hales et al, Phys. Rev. Lett. 85, 78 (2000); W. Ray 
et al, Phys. Rev. E 73, 026219 (2006). 

D. G. Luchinsky and P. V. E. McClintock, Nature 389, 
463 (1997). 

M. Morillo, J. M. Casado, and J. Gomez-Ordonez, Phys. 
Rev. E 55, 1521 (1997). 

M. I. Dykman et al, Phys. Rev. Lett. 68, 2718 (1992). 
D. Ludwig, SIAM Rev. 17, 605 (1975). 
H. B. Chan and C. Stambaugh, Phys. Rev. Lett. 99, 
060601 (2007). 

L. Onsager and S. Machlup, Phys. Rev. 91, 1505 (1953); 
S. Machlup and L. Onsager, Phys. Rev. 91, 1512 (1953). 



